arX1v:2109.07511v2 [astro-ph.HE] 28 Feb 2022 


Black hole spin-orbit misalignment in the x-ray binary 
MAXI J1820+070 


Juri Poutanen!**, Alexandra Veledina!*:?, Andrei V. Berdyugin!, 

Svetlana V. Berdyugina^, Helen Jermak?, Peter G. Jonker®’, Jari J. E. Kajaval^?, 

Ilia A. Kosenkov!, Vadim Kravtsov!, Vilppu Piirola!, Manisha Shrestha??, 
Manuel A. P. Torres!?:!!, Sergey S. Tsygankov!? 


1: Department of Physics and Astronomy, FI-20014 University of Turku, Finland 

2: Space Research Institute (IKT) of the Russian Academy of Sciences, 117997 Moscow, Russia 
3: Nordic Institute for Theoretical Physics (Nordita), KTH Royal Institute of Technology and 
Stockholm University, SE-10691 Stockholm, Sweden 

4: Leibniz-Institut für Sonnenphysik, 79104 Freiburg, Germany 

5: Astrophysics Research Institute, Liverpool John Moores University, L3 5RF Liverpool, 
United Kingdom 

6: Department of Astrophysics, Institute for Mathematics, Astrophysics and Particle Physics 
(IMAPP), Radboud University, NL-6500 GL Nijmegen, The Netherlands 

7: Space Research Organisation of the Netherlands (SRON), Netherlands Institute for Space 
Research, NL-2333, CA Leiden, The Netherlands 

8: Centro de Astrobiología, Villanueva de la Cafiada, S-28692 Madrid, Spain 

9: Department of Physics & Astronomy, University of Denver, Denver, CO 80208, USA 

10: Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain 

11: Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, 
Spain 


*Corresponding author. E-mail: juri.poutanen Q utu.fi. 


The observational signatures of black holes in x-ray binary systems depend 
on their masses, spins, accretion rate and the misalignment angle between the 
black hole spin and the orbital angular momentum. We present optical po- 
larimetric observations of the black hole x-ray binary MAXI J1820+070, from 
which we constrain the position angle of the binary orbital axis. Combining 
this with previous determinations of the relativistic jet orientation axis, which 
traces the black hole spin, and the inclination of the orbit, we determine a 
lower limit of 40? on the spin-orbit misalignment angle. The misalignment 
has to originate from either the binary or black hole formation stage. If other 
x-ray binaries have similarly large misalignments, these would bias measure- 
ments of black hole masses and spins from x-ray observations. 
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Black holes can be characterized with just two parameters: mass and spin. When a black 
hole resides in a binary system, accreting material from a companion donor star through the 
accretion disk, there are additional parameters that determine its observational signatures: the 
mass accretion rate and the misalignment angle between the black hole spin and the orbital 
axis. Standard methods to measure black hole spin from x-ray observations — iron line spec- 
troscopy (/) or modeling of the accretion disk spectrum (2) — assume that the misalignment 
angle is small. Conversely, the standard interpretation of low-frequency quasi-periodic oscilla- 
tions in x-ray and optical observations of black hole x-ray binaries as precession of the accretion 
disk (3—5), requires the assumption that the misalignment angle is non-negligible. Substantial 
misalignment is theoretically predicted for x-ray binaries that received high velocities during 
formation (6). The misalignment angle must be inherited from the formation process, because 
it can only decrease when the black hole is accreting (7). Gravitational wave observations of 
merging black holes have detected signatures of orbital precession (8) indicating non-zero mis- 
alignment in these systems (9), though they might not be representative of the wider population. 

Measuring the misalignment angle in x-ray binaries requires determining the three-dimensional 
orientation of the black hole spin and orbital axis. Accreting black holes often show relativistic 
jets, which are launched along an axis determined by the black hole spin direction (70). The jet 
inclination angle can be directly obtained in some cases from radio observations (//), whereas 
the jet position angle can be measured using either radio or x-ray imaging. Combining these 
two angles allows the black hole spin orientation to be determined. Orbital parameters, such as 
period and orbital inclination, can be determined using spectroscopic measurements of radial 
velocities of the donor star taken during quiescence, the stage at which accretion to the black 
hole is reduced and optical emission is not dominated by the accretion disk, through orbital 
modulation of the optical photometry, and using constraints from the presence or absence of 
x-ray and optical disk eclipses (72). 

The black hole x-ray binary MAXI J1820+070 was discovered as a transient x-ray source 
on 2018 March 11 (13). X-ray quasi-periodic oscillations, detected shortly after the discov- 
ery, were observed for more than 100 days (74). Ejections of material traveling at relativistic 
velocities have been observed from this source in both radio and x-rays, indicating the jet in- 
clination (measured from the line-of-sight) is tje = 63° + 3° and the position angle (measured 
on the plane of the sky from North to East) is Oje = 25°1 + 174 (15-17). Both angles were 
determined to be stable over the observed duration of the outburst. The orbital inclination has 
been constrained to the range 66° < torp < 81° by the lack of x-ray eclipses and the detection 
of grazing optical eclipses (72). To determine the orientation of the orbital axis requires one 
further parameter: the orbital position angle bor». 

We monitored MAXI J1820+070 in the optical B, V, and R photometric bands using the 
Double Image Polarimeters (78, /9) during the 2018 outburst and in quiescence. We obtained 
the source intrinsic linear polarization by subtracting the foreground interstellar polarization, 
measured from nearby field stars. During the outburst, when the relativistic jets were detected 
at radio frequencies, the intrinsic linear polarization degree (PD) in the V and R bands reached 
0.5% at a polarization angle (PA, also measured from North to East) of 23^—24^, which co- 


incides with the jet position angle within the uncertainties (20, 27). After the source faded in 
the x-rays, the PD increased by a factor of 5-10 and the PA changed by 40°+4° to —17°+4° 
(Table S1 and Fig. 1) (23). This increase in PD is most prominent in the B-band, which also has 
the highest PD in the range 1.5-5%; the R-band polarization changes from 0.4% to 2%. The 
PA is most precisely determined in the B-band, which also shows the least variability, with the 
mean being (PA) = —19°7 + 1°2. 


PD (%) 


PA (deg) 


0.0 0.2 0.4 0.6 0.8 1.0 
Phase 


Fig. 1. Observed optical polarization properties of MAXI J1820+070. (A) Intrinsic polar- 
ization degree (PD) and (B) polarization angle (PA) of MAXI J1820+070 in quiescence are shown as a 
function of orbital phase [using a published ephemeris (22)]. The intrinsic values were obtained from the 
observed ones by subtracting the foreground interstellar polarization, measured from nearby field stars. 
Blue circles, green triangles and red squares correspond to B, V, and R bands, respectively, with the 
error bars showing the 68% confidence level. Polarization is strongest in the B band and weakest in the 
R band, while the angle does not change substantially. 


We identify three properties of the quiescent-state polarization: it is strongest in the blue 
part of the optical, with approximate dependence on frequency v as PD(v) œ v’ (Table S1, 
Fig. 2), the PD remains high in the range 0.5-5% and the PA is stable. The PA undergoes appar- 
ently stochastic variations with an amplitude of <10° with no dependence on the orbital phase. 
These properties constrain the mechanism of the polarized emission. We modeled broadband 
photometric data obtained with the Liverpool Telescope and the Swift Ultraviolet and Optical 
Telescope (UVOT) together with the polarized fluxes (Fig. 2). We decompose the total spectral 
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energy distribution into three components: a companion star (contributing ~ 25% to the R-band 
flux), an accretion disk with inner temperature Ta ~ 6, 200 K and inner radius R4 z 6x 10? cm, 
and an additional ultraviolet (UV) component with blackbody temperature Tjj, œ~ 15, 000 K and 
radius Ry, z 9 x 10? cm (Table S4). The properties of the polarized flux are consistent with 
being produced by the UV component with constant PD of 5-846. 
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Fig. 2. Spectral energy distribution (SED). The average SED of MAXI J1820+070 (red dia- 
monds) as observed by the Liverpool Telescope and Swift/UVOT telescope in July 2020 and corrected 
for reddening with color excess E(B — V) = 0.29. The photometric bands are indicated at the top 
of the figure. The black dotted lines give the lower/upper limits on the flux for lower/higher extinction 
with E(B — V) = 0.25 and 0.325, respectively. The polarized flux divided by the best-fitting model 
polarization degree Puy = 0.055 (i.e. multiplied by a factor of ~18) is shown with blue triangles. Error 
bars show 68% confidence level. The black line gives the total model flux consisting of the companion 
star modeled as a blackbody (pink dot-dashed line), accretion disk (red dotted) and a hot blackbody (blue 
dashed). The spectrum of a K7 star (24) is shown as the solid green line for comparison. 


The jet cannot be the source of the polarized emission because its optically thin synchrotron 
spectrum is red, inconsistent with the observed blue spectrum of polarized light. Moreover, 
the PA is offset by about 40° from the jet position angle. The absence of detectable orbital 
variations in the PA excludes a hot spot origin. An optically thick accretion disk is excluded by 
the high PD and blue spectrum. A potential source of the polarized emission is scattering of the 
accretion disk radiation in the hot optically thin and geometrically thick accretion flow close to 
the disk inner radius (23, 25), which can also be responsible for the observed UV excess. This 
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mechanism would produce polarization parallel to the meridional plane, i.e. the plane formed 
by the orbital axis and the direction towards the observer. Another possibility is dust scattering, 
thought to be responsible for blue polarized spectra observed from accretion disks around some 
supermassive black holes (26). The presence of dust in quiescent-state black hole x-ray binaries 
has been inferred from the detection of the mid-IR excess in two systems (27). If dust is located 
within a flattened envelope, in the wind around the accretion disk, or in a circumbinary disk, the 
resulting polarization vector would also be parallel to the meridional plane. However, if dust 
forms an extended, approximately spherical structure at high elevation above the accretion disk, 
the polarization would be perpendicular to the meridional plane. We consider the latter scenario 
to be implausible, as a nearly spherical envelope cannot produce the high observed PD. A dust 
scattering mechanism would not explain the UV excess, because the disk does not emit in that 
range and hence there are no photons to be scattered by the dust. 
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Fig. 3. Geometry of the system from the observer’s perspective. The grey plane is the plane of 
the sky, labelled with North and East axes, perpendicular to the line of sight towards the observer 6. The 
angles between the line of sight and the vectors of the orbital angular momentum €) and the black hole 
spin are the inclinations torp and ipn. The corresponding position angles 05,5 and Opn are the azimuthal 
angles projected onto the sky, measured from North to East. The misalignment angle f is defined as the 
angle between 8 and Ê. The red cone indicates the jet and the blue ellipse indicates the companion star 
orbit around the black hole, which is at the coordinate center. 


Independent of the spectral modeling and geometry of the emission, the stability of the 
PA (most evident in the 5-band, Fig. 1) over the orbital phase suggests that the polarization 
is related to the orbital axis, either parallel or perpendicular to it. Hence, the observed PA 
provides information about the position angle of the orbital axis. The misalignment angle 5 can 
be determined from 


cos B = COS ibh COS Zorh + SiN dy sin torb cos A, (1) 


where tpn is the inclination of the black hole spin vector (measured from the line-of-sight) and 
A = bn — lor is the difference between the position angles of the black hole spin vector Opn and 
the orbital angular momentum ĝo» (the geometry is illustrated in Fig. 3). If the black hole spin 
vector is directed along the southern approaching jet, then its inclination tpn = tjet = 63° + 3° 
and its position angle is #4, = 180° + 6j, = 205?1 + 1°4 (15-17). The smallest misalignment 
B =~ 42° is achieved when the orbital spin is also directed South at 05, = (PA)+180°=160°3 + 
172 (because the PA has an ambiguity of 180°), at the inclination torb œ~ 73°. The probability 
distribution for £ in this case is shown in Fig. 4. The radial velocity measurements (/2) do not 
differentiate between orbital inclinations torp and 180? — torn, so there is a second solution with 
ig ^£ 107? and 8 ~ 63°. If either the orbital angular momentum or the black hole spin is 
instead directed to the North, the black hole rotation is then retrograde, resulting in 6 ~ 117° 
or 138? for the same two solutions for the orbital inclination as above. 

If the polarization vector is perpendicular to the meridional plane, the orbital position angle 
can take values orp = (PA) -- 90? or (PA) + 270°. This geometrical arrangement leads to nearly 
identical values for 3 because the difference between jet position angle and (PA) is about 45°. 
All possible cases for the orientations of the black hole and orbital spins, the resulting values 
for f and the azimuthal angle of the black hole spin in the orbital plane are listed in Table S5. 
Corresponding probability distributions are shown in Figs. $4 and S5. 

The difference of ~ 45? between the jet position angle and the PA indicates = 40° misalign- 
ment between the black hole spin and the orbital angular momentum. This result is independent 
of modeling or geometric ambiguities, because it relies only on the observed difference between 
the polarization angle and jet position angle. 

During outbursts, when the matter reaches the black hole, this misalignment affects the 
innermost regions of the accretion disk. For a non-zero spin, particles moving around the black 
hole in orbits tilted with respect to the black hole equatorial plane undergo precession at a 
rate that decreases with radius (3). Hence, a tilted disk is subject to twist and warp. A high 
misalignment adds complications to the models of quasi-periodic oscillations observed in black 
hole x-ray binaries, which rely on precession of the inner parts of the accretion flow, implying 
the whole flow is misaligned by 20 from the orbital axis in some phases (3). For 5 ~ 40° the 
inner parts of the accretion disk would need to become almost perpendicular to its outer parts. 
Most models assume smaller misalignment angles, typically, 6 ~ 10° — 20° (3-5), although 
highly inclined possibilities with 6 ~ 45? — 65° have sometimes been considered (26). 
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Fig. 4. Probability distribution function for the misalignment angle. The distribution normal- 
ized to the peak value is shown for the smallest misalignment angle possible. This case corresponds to 
the black hole spin directed along the southern approaching jet and the orbital spin being directed South 
at a position angle 05,5, = (PA) + 180° and inclination forb zz 73°. The red hatched region corresponds 
to the 68% confidence interval (i.e. between 16th and 84th percentiles of the posterior probability distri- 
bution). Distributions of 8 for the other seven possible combinations of 05,5, torb and ipn are shown in 
Fig. S4. 


High misalignment has previously been suggested based on observations of the gamma- 
ray light curves produced by the jet in Cyg X-3 (29), and differences between orbital and jet 
inclination angles are 15° in GRO J1655-40 (7) and 50° in V4641 Sgr (30) though the latter 
is highly uncertain. Misalignment has also been theorized, based on the inferred high kick 
velocities of x-ray binaries acquired during formation (6). For the black hole x-ray binary 
MAXI J1820+070, the high misalignment was identified only after obtaining the constraints on 
the position angle of the orbital angular momentum 0,,,. Without information on the binary 
plane orientation, we would have obtained only a lower limit on the misalignment angle in 
MAXI J1820--070 of 25°, because the orbital inclination is only marginally different from the 
jet inclination. 

Our results demonstrate the need to treat the misalignment angle as a free parameter when 
measuring black hole masses and spins. Assuming that the black hole spin and the orbital 
angular momentum are aligned introduces a systematic bias on measurements (12, /5, 31). A 
large misalignment angle is expected to drive precession of the binary orbital plane, altering 
the gravitational waves emitted during a subsequent merger event (9). Evidence for orbital 
precession has been found from population properties of black hole mergers observed using 
gravitational waves (6). 
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Materials and Methods 


DIPol-UF/NOT observations 


We have performed polarimetric observations of the black hole x-ray binary MAXI J1820+070 
(13,34, 35) both during the bright stages of the outburst in 2018 as well as in the quiescence. The 
data during the outburst (20, 27) were obtained with the Double Image Polarimeter-2 (DIPol- 
2) polarimeter (78) mounted on the Tohoku 60-cm telescope (T60) at Haleakala observatory, 
Hawaii. During the quiescence MAXI J1820+070 was observed with the DIPol Ultra Fast 
(DIPol-UF) (19) at the 2.56-m Nordic Optical Telescope (NOT) in July 2019, April and July 
2020, and July 2021 (Table (32). Both polarimeters have similar design: incoming light 
passes through a modulator (superachromatic half-wavelength plate, HWP), it is then split by 
the plane-parallel calcite plate into two orthogonally polarized rays (ordinary and extraordinary) 
and further is reflected by two dichroic beam splitters to produce o- and e-images on three 
charge-coupled devices (CCDs) simultaneously in the B, V, R bands. This design optically 
eliminates the sky polarization at the registration stage (even if it is variable), reaching up to 
10^? polarization precision (/8, 19, 36). The instrumental polarization is below 1074 level and 
can be accounted for using zero-polarization standards (/9) but is negligible for our study. The 
zero point of the PA was determined by observing high polarization standards HD 161056 and 
HD 204827 (19). 

For each independent measurement of linear polarization, four sequential images with the 
HWP rotated by 2275 in each step are obtained. Each measurement cycle consists of 16 images 
recorded simultaneously by three CCDs for one full rotation of the HWP modulator (0?—360?) 
giving four independent measurements of the normalized Stokes parameters (q, u) in the D, V, 
and R bands. This algorithm helps eliminate effects arising from dust particles on the retarder, 
nonparallelism of rotating components, etc. The images produced by DIPol-UF require standard 
CCD calibrations (36). 

The average normalized Stokes parameters (q, u) are then obtained from 16 to 40 individual 
measurements using a 2c weighting algorithm (/9, 37). The typical duration of an observational 
window where the data were averaged was one hour. The polarization produced by the interstel- 
lar medium (ISM) has been estimated from numerous (more than 400) observations of sample 
of field stars (stars 42, 3, 6, 7, 9 from fig. 9 of ref. (20)), which are close in distance to the 
target as indicated by their parallaxes (38). These normalized Stokes parameters (qs, Uism) 
were subtracted from the measured values of the normalized Stokes parameters to obtain the 
intrinsic polarization information of the source (intr, Uintr). We then define the complex linear 
polarization quantity as P = dintr + 72Uintr. Intrinsic polarization degree (PD) P and polarization 
angle (PA) 0 are then obtained from the formulae 


la LE 1 
P — |P| m be + Ukat 0 = 2 arg(P). (S1) 


Because our PD measurements typically have significance 5—10c, the bias in the measured PD 
is negligible and the uncertainty in PD AP is the same as the uncertainty on individual Stokes 
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Figure S1: Variation of flux and polarization. (A) The dependence of the observed flux (not 
corrected for reddening), (B) intrinsic source PD and (C) PA on time during the July 2020 
observations of MAXI J1820+070 with DIPol-UF in three filters: B (blue circles), V (green 
triangles), and R (red squares). (D) Dependence of the PD on flux in the same three filters. An 
anticorrelation is evident. 


parameters. The uncertainty of the PA is estimated as A0 — AP/(2P) (39). The intrinsic PD 
and PA computed following this procedure are reported in Table [S1] 

Polarization observed during quiescence shows no clear dependence on the orbital phase 
(Fig. 2), apart from some spread of the points in close orbital phases. The PD has a blue spec- 
trum depending on frequency v as x v3, which is inconsistent with the red spectra expected 
from a jet or accretion disk. The PA is very stable in the B-band, where the PD is the highest 
and the relative uncertainty is smallest. Therefore, we computed the average PA in that band as a 
inverse-variance weighted mean of individual PAs. Since the deviations of individual measure- 
ments exceed the measurement uncertainties, the standard error of the weighted mean was cor- 
rected for overdispersion by the square root of the reduced x, resulting in (PA) = —19°7+1°2. 
We also computed the mean and its uncertainty using a bootstrap (sampling with replacement) 
method (40) arriving at an identical value. The mean PA in the V- and R-bands have similar 
values, but larger uncertainties (Table S1). 

Using the images from DIPol-UF we also measured photometry of MAXI J1820+070 rela- 
tive to two nearby stars (#1 and 2 from table 2 in ref. (20)). Absolute photometry during July 
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Table S1: Intrinsic polarization of MAXI J1820+070. The intrinsic values of the PD and PA 
derived from observations of MAXI J1820+070 during the outburst with DIPol-2 (20, 27) and 
during quiescence with DIPol-UF at NOT. Uncertainties are 1o. The Stokes parameters have 
been corrected for the interstellar polarization. The dates for observations during quiescence 
correspond to the middle of observing windows of a typical duration of one hour. 


Date B V R 


(MJD) PD (%) PA (deg) PD (%) PA (deg) PD (%) PA (deg) 


Outburst 
58195-58222 0.28 + 0.01 9.2 + 1.0 0.36 + 0.01 22.9+ 1.0 0.30 + 0.01 29.0 + 0.9 
58223-58234 0.34 + 0.02 8.8 1.4 0.51 + 0.02 23.4+ 1.4 0.53 + 0.02 23.9+1.1 
58312-58344 0.16 + 0.01 15.8 + 1.6 0.15 + 0.01 13.4+ 2.3 0.02+0.01 39.1+11.0 
58406-58428 0.06 + 0.04 3.0 + 15.4 0.13 + 0.06 2.8+ 12.4 0.09+ 0.04 44.8 + 12.5 
Quiescence 
58686.93 3.86 + 0.29 23.8 + 2.1 1.81 +0.24 —23.8 + 3.8 0.70 + 0.14 24.1 4: 5.7 
58688.01 2.35 + 0.38 25.9 + 4.7 1.09 + 0.25 8.0 € 6.5 1.70 + 0.10 42.4 X 1.7 
58961.15 1.47 + 0.28 26.8 + 5.4 0.52 +0.24 —4.2 + 13.3 0.63 + 0.19 8.5+ 8.7 
58964.16 4.31 + 0.28 18.5 + 1.8 3.10 +0.29 —17.5 € 2.7 1.05 «& 0.18 —2.0 + 4.8 
58965.11 2.81 + 0.24 23.5 + 2.5 2.12+0.21 —35.1 + 2.9 0.84 + 0.18 26.3 + 6.2 
58967.12 4.53 + 0.28 24.1 1.8 2.09 +0.39 —11.0 5.3 1.09 + 0.31 27.0 + 8.3 
59050.94 4.58 + 0.34 15.3 + 2.1 3.57 +0.31  —16.5 + 2.5 1.54 + 0.18 14.5 + 3.4 
59050.98 5.03 +0.30 —11.0 1.7 3.38 +0.34  —11.7 4 2.9 1.62 + 0.16 0.34 2.9 
59051.05 3.35 +0.26 —23.5 + 2.2 2.04 + 0.23 —17.1 + 3.2 1.08 + 0.18 12.1 + 4.9 
59051.10 3.14 + 0.30 10.7 + 2.7 1.95 + 0.29 —4.8 + 4.3 0.95 + 0.16 5.6 + 4.9 
59051.92 4.08 + 0.27 18.1 + 1.9 3.40 + 0.25  —12.9 + 2.1 1.39 + 0.15 18.0 3: 3.1 
59051.97 4.67 + 0.30 16.5 + 1.8 3.27 +0.30  —17.5 + 2.6 1.96 +0.15 —17.9 + 2.2 
59052.03 4.83 + 0.31 16.7 + 1.9 3.37 +0.24 —10.3 + 2.1 1.68 + 0.15 20.0 + 2.6 
59052.09 3.78 + 0.23 21.5+1.7 2.47 +0.22  —28.0 + 2.5 1.06 + 0.16 12.4 + 4.3 
59052.92 4.14 + 0.23 24.2 + 1.6 2.46 + 0.23  —16.6 + 2.7 1.21 + 0.14 18.4 + 3.3 
59052.98 3.30 + 0.28 20.5 + 2.4 2.18 +0.23  —23.6 + 3.0 1.18 +0.17 —14.0 + 4.0 
59053.92 4.02 + 0.25 12.8 + 1.8 2.22 + 0.26 —3.5 + 3.4 1.24 + 0.15 0.2 + 3.4 
59053.97 4.30 + 0.27 16.5 + 1.8 3.18 +0.24  —10.4 2.1 1.51 +0.15 —13.1 + 2.9 
59400.99 1.62 + 0.12 30.5 + 2.0 0.77 +0.10 —17.1 3.9 0.41 0.07 6.4 3: 5.1 
59401.94 2.86 + 0.21 25.4 + 2.1 2.04+0.20  —29.9 + 2.8 0.36+0.14 —6.1 + 10.7 
59051—59054 4.00 + 0.19 17.5 + 1.3 2.65 +0.19  —14.7 € 1.9 1.32 +0.08 —12.8 + 2.3 
58686-59402 3.18 + 0.22 19.7 + 1.2 1.89+0.19 —16.8 + 1.9 0.94 + 0.11 18.6 + 3.3 
Interstellar polarization 
58195-59054 0.81 + 0.03 64.0 + 1.1 0.71 + 0.03 68.4 + 1.2 0.60 + 0.02 64.4 + 0.8 
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Table S2: Observed flux of MAXI J1820+070. Mean fluxes and their 1c uncertainties derived 
from observations of MAXI J1820+070 during quiescence in July 2020 with DIPol-UF at NOT. 


Date Flux (mJy) 

(MJD) B V R 
59050.94 0.0732 + 0.0014 0.1391 + 0.0019 0.1858 + 0.0020 
59050.98 0.0747 + 0.0022 0.1495 + 0.0030 0.1966 + 0.0029 
59051.05 0.1000 + 0.0022 0.2002 + 0.0034 0.2491 + 0.0035 
59051.10 0.1090 + 0.0018 0.2142 + 0.0030 0.2772 + 0.0028 
59051.92 0.0819 + 0.0015 0.1577 + 0.0021 0.2082 + 0.0023 
59051.97 0.0700 Æ 0.0014 0.1420 + 0.0022 0.1901 0.0022 
59052.03 0.0681 + 0.0015 0.1395 + 0.0024 0.1941 0.0025 
59052.09 0.0813 + 0.0018 0.1572 + 0.0026 0.2142 + 0.0027 
59052.92 0.0902 + 0.0015 0.1745 + 0.0023 0.2255 + 0.0023 
59052.98 0.0922 + 0.0021 0.1754 0.0031 0.2323 + 0.0031 
59053.92 0.0861 + 0.0017 0.1759 + 0.0027 0.2275 + 0.0028 
59053.97 0.0736 + 0.0013 0.1512 + 0.0020 0.1995 + 0.0020 


2020 observations was obtained using simultaneous observations by the Liverpool Telescope 
(LT), see Fig. and Table [52| The flux of MAXI J1820+070 in all bands is anti-correlated 
with the PD (Fig. (SID). This anti-correlation together with the spectral dependence of the PD 
can be interpreted as an interplay between two components: one, polarized, with contribution 
growing towards blue and another, unpolarized or, possibly, weakly polarized, variable com- 
ponent dominating the flux in the red. Following the spectral decomposition (Fig. 2), the UV 
component produces only 20% of the total flux in the R band. This makes PA measurements of 
the UV component less reliable in this band. Stronger variations of the PA in may be caused 
by the presence of a red and polarized at a level of 0.1—0.396 component, whose PA is different 
from that of the UV component. Such red and polarized component, with the PA consistent 
with the jet position angle, has been observed during the outburst (Table|S1). 


Liverpool Telescope observations 


Photometric observations in the optical band were performed using the optical imaging com- 
ponent of the Infrared-Optical suite of instruments (IO:O) with the Sloan Digital Sky Survey 
(SDSS) i'r'z', Bessel V and B filters on the LT (4/7) on La Palma, Spain. The object was 
observed during four nights from 2020 July 20 to 24, i.e. Modified Julian Dates (MJD) 59050- 
59054. LT image reduction was provided by the basic IO:O pipeline (42). As the photometric 
standards, we used four stars with known Panoramic Survey Telescope and Rapid Response 
System (PS1) point spread function magnitudes. To convert PS1 magnitudes to SDSS and 
Johnson-Cousins magnitudes we used equation (6) and table 6 from ref. (43). To obtain fluxes 
of the object from its magnitudes we used standard zero-points for SDSS and Johnson-Cousins 
systems (44, 45). The fluxes and their uncertainties are presented in Table [S3] 
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Table S3: Log of observations with the Liverpool Telescope. The observed fluxes of 
MAXI J1820+070 and their 1c uncertainties in various filters. 


Date Filter Exposure Flux Date Filter Exposure Flux 

(MJD) (s) (mJy) (MJD) (s) (mJy) 
59050.92075 i 100 0.2700+0.0014 59052.92842 z 100 0.3944+0.0033 
59050.92213 a 100 0.2468+0.0015 59052.92980 z 100 0.40370.0033 
59050.92350 i 100 0.2482+0.0015 59052.93193 V 200 0.1605+0.0011 
59050.92505 z 100 0.3990+0.0028 59052.93463 B 200 0.09740.0007 
59050.92642 z 100 0.4093+0.0029 59052.93637 r 40 0.3063+0.0018 
59050.92780 z 100 0.3733+0.0029 59052.93705 T 40 0.2142+0.0018 
59050.93052 V 300  0.14402-0.0009 59053.92120 a 100  0.3071-:0.0017 
59051.01279 i 100 0.3321+0.0014 59053.92257 D 100 0.3600+0.0018 
59051.01417 D 100 0.3075+0.0014 59053.92395 a 100 0.3450+0.0018 
5905 1.01554 i 100 0.3164+0.0014 59053.92550 z 100 0.39170.0033 
59051.01709 z 100 0.42080.0027 59053.92687 z 100 0.3598+0.0033 
59051.01846 A 100 0.50480.0028 59053.92825 z 100 0.3736+0.0034 
5905 1.01984 z 100 0.4703+0.0027 59053.93038 V 200 0.1369+0.0011 
59051.02254 B 300  0.1263::0.0005 59053.93308 B 200 0.0762+0.0007 
59051.09342 D 100 0.3476+0.0016 59053.93482 r 40 0.2281+0.0019 
5905 1.09479 i 100 0.3489+0.0016 59053.93550 r 40 0.2174+0.0019 
59051.09617 D 100 0.3809+0.0015 59053.95474 a 100  0.2752::0.0016 
59051.09771 z 100 0.4877+0.0032 59053.95612 a 100  0.2927-:0.0016 
59051.09909 z 100 0.49630.0033 59053.95750 a 100 0.2633+0.0017 
5905 1.10047 z 100 0.4706+0.0033 59053.95904 z 100 0.3587+0.0032 
59051.10168 r 40 0.2339+0.0019 59053.96042 z 100 0.3824+0.0030 
59051.10236 r 40 0.2592+0.0018 59053.96180 z 100 0.40970.0033 
59051.92206 D 100 0.3003+0.0015 59053.96393 V 200 0.1467+0.0010 
59051.92343 D 100 0.2982+0.0015 59053.96663 B 200  0.07993-0.0007 
59051.92481 i 100 0.2860+0.0015 59053.97116 a 100 0.3032+0.0042 
59052.01060 D 100 0.2838+0.0016 59053.97254 a 100 0.3390+0.0016 
59052.01197 i 100 0.2693+0.0015 59053.9739 1 a 100 0.3308+0.0016 
59052.01335 D 100 0.2501+0.0016 59053.97546 z 100 0.4278+0.0032 
59052.09336 D 100 0.2979+0.0015 59053.97684 z 100 0.40170.0032 
59052.09474 i 100 0.2800+0.0016 59053.97821 z 100 0.4193+0.0031 
59052.0961 1 i 100 0.2974+0.0015 59053.98035 V 200  0.15352-0.0010 
59052.92275 i 100 0.3161+0.0017 59053.98305 B 200 0.0601+0.0007 
59052.92412 D 100 0.3058+0.0016 59053.98479 f 40 0.2427+0.0017 
59052.92550 D 100 0.3048+0.0017 59053.98546 r 40 0.2100+0.0017 
59052.92705 z 100 0.3686+0.0034 
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Swift/UVOT observations 


The Neil Gehrels Swift Observatory (Swift) (46) observed MAXI J1820+070 in the low state 
with the X-ray Telescope (XRT) and Ultraviolet/Optical Telescope (UVOT) instruments (partly 
quasi-simultaneously with DIPol-UF/NOT) between 2020 July 20 — September 6. All the XRT 
data were taken in the photon counting mode with about 9 ks of total exposure. An averaged 
spectrum was extracted using the online tools (47) provided by the UK Swift Science Data 
Centre. The spectrum was rebinned to have at least 1 count per energy channel in order to apply 
W -statistic (45, 49). The data were fitted with a model consisting of a power law modified 
by the interstellar absorption TBABS*POWERLAW in the 0.5-10 keV band using the XSPEC 
package (50). Fixing the hydrogen column density to Ny = 1.6 x 10?! cm? (57) we obtained 
a power-law photon index I' = 1.62 + 0.24 with W -statistics of 42.7 for 46 degrees of freedom. 
We find flux corrected for absorption of (3.3 + 0.7) x 10713 ergcm ?s^! in the 0.5-10 keV 
energy range. 

The UVOT data were processed with the UVOTPRODUCT tool (52—54). We used a 5" source 
aperture around the known position of MAXI J1820+070 (38), and a circular source-free 10" 
background aperture about 40" east of it. After comparing the magnitudes obtained between 
July 20 — September 6 to other quiescence data taken when the x-ray counting rate was low, we 
found that during this period MAXI J1820+070 was roughly 0.5-1 magnitude fainter in all the 
UV filters than in the data taken in 2019 or earlier in 2020. The same trend was observed in the 
V -band and, to a lesser extent, in the R-band. Therefore, we decided to stack all the 2020 July 
20-September 6 UVOT data together to obtain time averaged fluxes in the u, uvw/, uvm2 and 
uvw2 filters for the spectral energy distribution (SED). Conversion of the background corrected 
count rates obtained from the UVOTPRODUCT to the fluxes is non-trivial, as it depends on the 
assumed spectral shape (53). Therefore, for reliable estimate of the fluxes we performed direct 
spectral fitting with a broad-band SED model using the XSPEC package (Fig. 2). 


Decomposition of the broadband spectrum 


The broad-band average spectrum obtained with the LT and Swift/UVOT is shown in Fig. 2. 
The spectrum was corrected for reddening using model of ref. (55) with extinction Ay = 0.89 
(corresponding to the color excess E(B — V) — 0.29) obtained from the hydrogen column 
density Ny = (0.16 + 0.02) x 10?2cm~? (51) using transformation Ay = 5.59 Ny/10?? cm~? 
(56). The optical/infrared SED obtained with the LT is very red and can be described by a joined 
contribution of companion star and multicolor accretion disk, but the UVOT data show a UV 
excess. The excess is also seen in the non-dereddened data and is not an artefact of dereddening 
procedure. A similar excess has been detected previously in quiescent-state observations of the 
black hole x-ray binary A 0620—00 (57), but its nature remained a mystery. 

The blue PD spectrum differs from the red (even after correcting for the reddening) spectrum 
of the infrared to UV continuum. This rules out accretion disk as the source of polarized light. 
We infer that the component producing the UV excess is responsible for polarization. To extract 
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the average shape of the polarized component, we calculated the polarized fluxes as 
(PE = Pike (S2) 


where k is the index corresponding to one of the BV R-bands, Pj is the average PD and Fk 
is the average flux in this band. If PD of the polarized component is wavelength-independent, 
then the computed PF directly replicates its total spectral shape, with the normalization being 
smaller by the polarization fraction. The average PD was computed from the average Stokes 
parameters (dint; Uintr) obtained in 12 individual measurements in July 2020 during quiescence 
with DIPol-UF (Table|S1). Because absolute photometry with DIPol-UF is not very reliable, the 
average flux was obtained from the LT observations. The averages were computed as a inverse- 
variance weighted mean of individual measurements and their standard errors were corrected 
for overdispersion by the square root of the reduced y?. 

Published spectroscopy obtained during quiescence indicates a K-type spectral type for the 
companion star, which contributed about 20% to the flux in the r-band (/2, 22). The average 
observed r-band flux during the spectroscopy was about 0.33 mJy, higher than the average flux 
of 0.24 mJy in our LT observations. Thus, the contribution of the companion star to the r- 
band flux rose to 27%. The radius of the companion star, which fills its Roche lobe, can be 
established from the measured orbital period P5, = 0.68549 + 0.00001 d, projected rotational 
velocity vrot sin tory = 84 + 5 km s^ 1, and orbital inclination i44 = 73? + 6? (72) 


Por Tro l Lor 
R, = Perot Sinton) _ (1 19 1008) R; = (8.3 + 0.6) x 10! om. (S3) 
27 Sin torb 


A moderately evolved late-type star of spectral class K7 can satisfy the constraint on radius 
and contribute about 2796 to the r-band flux. A star with higher temperature overpredicts the 
companion star contribution in the r-band. For the fitting purposes, we approximate the stellar 
spectrum (24) by the blackbody with T, = 4000 K and R, = R. = 8.3 x 101° cm (see Fig. 2). 

Using XSPEC v.12.11 (50), we modelled the total SED from LT and Swift/UVOT and, 
jointly, the three points of the PF. The total spectrum is described by the model REDDEN 
(BBODYRAD,+DISKBB+BBODYRAD,). The REDDEN model describes interstellar extinction 
(55). Component BBODYRAD, corresponds to the spectrum of the stellar companion, modeled 
as a blackbody of fixed temperature 7; and radius R,. The DISKBB corresponds to the mul- 
ticolor accretion disk. It has two free parameters: normalization, which is related to the inner 
radius Rin, and temperature at that radius Tin. Component BBODYRAD, corresponds to the UV 
excess that is modeled by a blackbody with two free parameters: temperature 7}, and radius 
Ry». The PF is modeled by the reddened second blackbody multiplied by the polarization 
fraction Pyy. Blackbody and disk normalizations were converted to radii assuming the dis- 
tance to the source of 2.96 kpc as determined by the radio parallax (/5) and disk inclination 
torb = 73° (12). 

The best-fitting model parameters with the corresponding uncertainties are listed in Ta- 
ble [S4] We find the accretion disk temperature of about 6200 K, which is very close to that 
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Table S4: Best-fitting parameters of the SED model. Distance of 2.96 kpc (15) and inclina- 
tion 45, = 73° (12) are assumed. The temperature and the radius of the blackbody approximat- 
ing the SED of the stellar companion were fixed. 


Parameter Value Units 
Ti 6200 15 K 
Rin (5.6227) x 101° cm 
Tis 14900+2300 K 
Rib (8.9423) x 10? cm 
T: 4000 K 
Ry 8.3 x 101° cm 
Poy 0.0557 5.011 


expected for the disk in quiescence (58). The inner radius of the disk Rin z 6 x 10!? cm, 
which is 30% smaller than the estimate of the circularization radius for the measured P,,1, and 
companion-to-black hole mass ratio q = 0.072 + 0.012 (72). The outer radius has to be at 
least 2-3 times larger to avoid underpredicting the red part of the spectrum. This is consis- 
tent with the expectation that the accretion disk size does not exceed the tidal radius of about 
2.4 x 10!! cm (58). 

The temperature of the additional blackbody is 7p» ~ 15, 000 K and the characteristic size 
Ry, ~ 9x10? cm. This reproduces the shape of the polarized flux SED and the UV excess. The 
PD of this component is Pyy = 0.055* 0077, corresponding to 5-8% intrinsic PD (on average) 
of the polarized component. The broad-band SED and the best-fitting spectral components are 
shown in Fig. 2. The computed values of PF were divided by the constant Pyy to demonstrate 
the fit with the UV excess component. 

For lower and higher values of the color excess of 0.25 and 0.325, corresponding to the lo 
uncertainties in Ny, the best-fitting spectral parameters change. The lower E(B — V) leads to 
the 8% smaller disk temperature Tin and 13% larger inner radius Rin, 8% smaller blackbody 
temperature Tp» and 3% larger radius Ry. For the higher E(B — V), the effect is opposite: Tin 
is 7% larger, Rin is 9% smaller, Tb» is 7% larger and Rp» is smaller by 3%. This uncertainty 
does not affect any of our conclusions. 


The source of polarized light and the nature of the spectral components 


In this section we seek to answer two questions: what is the nature of the UV excess and 
what is the source of polarized radiation? The accretion disk itself cannot be the source of the 
UV emission, because in quiescence its temperature is lower as we see from the red optical 
spectrum. Another possibility is a hotspot (or hot line), the place where the accretion stream 
hits the disk. This component could be responsible for the UV excess. However, whether it 
can also be the source of polarization is questionable. The high temperature of this component 
implies that the matter is ionized. The PD from the optically thick electron-scattering dominated 
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atmosphere (59, 60) depends on the cosine of the viewing angle u = cos torb approximately 
as 11.7% x (1 — u)/(1 + 3.5824) (61). For inclination angles permitted by the absence of 
eclipses (12), 454 < 81°, the expected PD is < 6%. A high implied PD, Pyy 225-896, is barely 
consistent with that. 

On the other hand, high PD can be produced by synchrotron radiation in the ordered mag- 
netic field. The blue spectrum indicates that we might see the optically thick part of this ra- 
diation, with the possible transition to optically thin part at v > 10° Hz. Self-absorption 
becomes important at the turn-over frequency (62-64), v, + 3 x 105 Bet (ry? sn, where Be 
is the magnetic field in units of 10° G, 7, is the Thomson optical depth of electrons emitting 
at the turn-over frequency, ^; is their Lorentz factor and p is the index of the power-law distri- 
bution of electron number density on the Lorentz factor y, dn./dy x y’. For the observed 
v ~ 10!° Hz, this requires highly opaque source with 7,57 ~ 1 and/or high magnetic field with 
Bg ~ 1, which would be expected in the bright hard state, but inconsistent with the relevant val- 
ues for the source in (near-)quiescence, when both optical depth and the magnetic field drop by 
2 to 4 orders of magnitude (65). The highest theoretically possible PD below the self-absorption 
frequency (66) of Pnax = 3/(6p + 13) ~ 10% requires highly ordered — and constant — mag- 
netic field during two years, July 2019—July 2021. Hence, we find the synchrotron radiation to 
be an implausible source of the observed polarized flux. 

Substantial polarization could instead be produced by electron scattering of radiation in an 
optically thin slab if the seed photons are injected along the slab plane (67). The PD— (1 — 
ui”) /(3 — u?) reaches 33% edge-on and is a weak function of orbital inclination (for iorb = 66°, 
ie.  < 0.41, the PD is larger than ~30%). Scattering in the Thomson regime means that 
the scattered radiation does not gain a systematic shift in energy with respect to the incident 
continuum, and hence the peak of the spectrum of the polarized component directly probes 
the characteristic energy and spectral shape of the incident radiation. Hence, the source of 
the incident light should have narrow, blackbody-like shape, with characteristic temperature 
T ~ 15,000 K. Because this temperature is much higher than the disk inner temperature, and 
the polarization angle is independent of orbital phase, the location and physical properties of 
such a component are unclear. 

The additional polarized component may instead arise from Compton up-scattering of soft 
photons by hot electrons, which is related to a systematic shift in energy. A hot, geometrically 
thick accretion flow can exist in low-luminosity black hole binaries and active galactic nuclei 
(25,65). Alternatively, a hot layer (corona) atop of the cold disk, heated by the unwinding stream 
of matter, could be present. In both cases, the hot medium is expected to be optically thin, and 
may Compton up-scatter the photons coming from the cold accretion disk. We calculate the 
spectra and polarization properties of Compton scattering using iterative scattering method (68). 
We assume the spectrum of incident, non-polarized photons is a blackbody with temperature 
equal to the disk inner temperature Tin = 6200 K. For these seed photons, to explain the peak 
of the first scattering order at ~ 10!? Hz (as demanded by the spectral decomposition), we need 
to have electrons with temperature of about 200 keV. For such a high temperature, the PD is 
expected to drop with respect to the maximal, Thomson-regime PD by a factor of 3 (69). In 
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Figure S2: Polarization properties of the hot accretion flow model. (A) Flux from the spec- 
tral components: blackbody (Tin = 6200 K, black solid line) and first Compton scattering of 
these photons by hot medium (red dashed line) with kT, = 200 keV and scale-height given by 
parameter cos Qi, = 0.5, and spectrum of the polarized flux of the first scattering component 
(blue dot-dashed line) that has been scaled up to match the low-energy part of the intensity spec- 
trum (divided by the PD at low energies). This shows that the polarized flux has the same shape 
as the scattered component flux, justifying our joint model fitting of the UV flux together with 
the polarized flux (shown in Fig. 2). Only results for inclination of torp = 73° are shown, as the 
dependence on the inclination in its allowed range is minor. (B) PD of the scattered component 
as a function of frequency is shown for different inclinations torp = 81° (blue dashed line), 73? 
(green dotted) and 66? (red solid). PD for energies below the peak of the scattered component 
is nearly frequency-independent. The observed strong dependence of the PD on frequency is 
caused by the dilution of the polarized scattered radiation by the unpolarized accretion disk. 
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the case of inner hot flow, only the photons travelling at inclinations a to the disk normal that 
are larger than some critical value a,,i, can interact with the hot matter. This limiting angle 
is related to the scale-height H of the hot medium cot o4, ~ (H/Rin). Increasing the flow 
scale-height leads to a drop of PD. In Fig. we show the resulting flux, polarized flux and 
PD spectra at different viewing angles 7,,, = 66°, 73° and 81°. In this case the polarization is 
parallel to the disk normal. We find that the PD of scattered radiation is nearly independent of 
the photon frequency up to its peak at v ~ 10? Hz and therefore the spectral shape of the flux 
and PF of the scattered component are nearly identical in the optical band. 

For the slab-corona geometry, an appropriate parameter determining PD is the Thomson 
optical depth of the hot slab tp. Only for rr = 0.5 a substantial polarization exceeding 5% 
can be produced and the dominant electric-field oscillations are perpendicular to the projection 
of the disk normal. However, Comptonization in a slab of rr = 0.1 and electron temperature 
kT, = 200 keV overproduces the observed x-ray flux. This would imply that a more probable 
model that is more consistent with the observed spectral energy distribution and polarization 
properties is scattering of the disk radiation in the inner hot flow. Such a geometry favors the 
measured PA being parallel to the orbital axis. 

Finally, the blue PD spectrum could be produced by dust scattering of the accretion disk 
radiation. This process is thought to be responsible for similarly blue polarization spectra 
observed from supermassive black holes in Seyfert galaxies and quasars (26, 70-74). The 
dust would likely be located in a flattened envelope (equatorial wedge) around the accretion 
disk (73), where it is shielded from the inner disk radiation, or in a circumbinary disk (27). In 
this case, the polarization vector lies in the meridional plane. If instead the dust has a more 
spherical distribution, the PD is expected to be smaller and the polarization to be perpendicular 
to the meridional plane. The dust scattering model, however, does not explain the UV excess. 


Geometry 


Here we define the coordinate systems and derive the formulae to compute the misalignment 
angle between the black hole spin and the orbital angular momentum, and the azimuthal angle 
of the black hole spin projection in the orbital plane. We consider a Cartesian system with 
the x — y plane coinciding with the orbital plane. Thus the unit vector of the orbital angular 
momentum is Ê = (0,0, 1). We choose the direction to an observer to lie in the x — z plane 
at inclination angle torb as measured from the orbital axis (see Fig. [53] for geometry), so the 
Observer unit vector is 

6 = (sin torb, 0, cos Zorp). (S4) 


We assume that the black hole spin is directed at an angle 6 from the z-axis at azimuthal angle 
bn, which is measured from the x-axis in the counter-clockwise direction in the x — y plane as 
viewed from the top. The unit vector of the black hole spin is 


§ = (sin B cos Oy, sin B sin y, cos £). (S5) 
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Figure S3: Geometry of the system from the point of view of the binary. The orbital plane 
here coincides with the x — y plane of the Cartesian system with the orbital axis being along the 
z-axis. The observer is situated in the x — z plane at inclination angle torp as measured from the 
orbital axis. The black hole spin is directed at angle 6 from the z-axis at azimuthal angle ®,), as 
measured from the x-axis in the counter-clockwise direction in the x — y plane as viewed from 
the top. The red cones indicate the jet, and the blue disk indicates the binary orbit. Fig. 3 shows 
the same geometry from the observer's perspective. 
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We define the polarization basis with the unit vector €; directed on the sky in the direction 
of the projection of the orbital spin, vector és being directed to the left on the sky, and the third 
vector coincides with the observer direction: 


Q — cos iq Ô 


ĉi = s = ( COS torb, 0, sin tors); 
orb 
i (S6) 
e2 = (0, E 0), 
Êz = ô 
In this vector basis, the black hole and the orbital spin vectors can be represented as 
& = (sinipn cos A, sin ip, sin A, cos iph), (S7) 
Q = (siniorb, 0, cos torb), (S8) 


where A = bn — Dorp is the difference in position angles of the black hole and the orbital spins. 
The misalignment angle 8 between the black hole and the orbital axis is then given by the scalar 
product: 

cos B = Ô- ê = cos ibh COS orb + SIN iph SIN torb cos A. (S9) 


The direction cosines of the black hole spin vector § in polarization basis are 


$-6; = sin ipn cos A = cos D sin torb — sin B cos torb cos ph, 
$-65 = sin ipp sin A = — sin B sin Oy, (S10) 
& - êg = cos iph = cos D cos tory + sin B sin torb cos Ppp, 


allowing us to obtain the azimuthal angle of the black hole spin: 


Sin torb COS ibh — COS torb SIN iph COS A 
sin f É 
sin ipn sin A 
sin 8 
We assume that the black hole spin is aligned with the jet. If the spin is directed along the 


southern approaching jet, then its inclination ijn = tje = 63° + 3° and its position angle is 
Oban = 180° + 0j. = 20521 + 1°4 (15-17). If, on the other hand, the black hole spin points 


COS Oy. = 
(S11) 


sin Oy. = 


along the northern receding jet, then ibn = 180° — ije = 117° + 3° and its position angle is 
Obh = ig = 2521 E 124. 
The polarimetric data provide us with the average polarization angle (PA) = —19°7 + 


1°2, which carries information about orientation of the orbital axis on the sky. If dominant 
oscillations of the electric field lie in the meridional plane formed by the orbital spin and photon 
propagation direction, the position angle of the orbital spin can be either 05, = (PA) or (PA) + 
180°. The electric field oscillations can also be perpendicular to the meridional plane, then 
the orbital spin position angle is (PA) + 90° or (PA) + 270°. Furthermore, the radial velocity 
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Table S5: Geometrical parameters. Geometrical parameters for the 16 possible cases identi- 
fied by letters A—P of relative orientation of the orbital and the black hole spins. 


ibh ijet = 63° 3° 180° — ijet = 117° 3° 
Obh Pet + 180° = 205° 1 + 1°4 Ojo = 25° 1 + 174 
torb 73° + 6° 107° + 6° 73° + 6° 107° + 6° 
Polarization parallel to the meridional plane 
borb = (PA) B (deg) 117.34 4.3 B 137.6 + 2.4 C 62.7 + 4.3 D 42.4 + 2.4 
=—19°7 + 1?2 Phn (deg) 45.2 + 3.4 69.7 + 7.0 225.2 + 3.4 249.0 + 7.0 
= 5 B (deg) 42.4 + 2.4 62.7 + 4.3 137.6 + 2.4 117.34 4.3 
Pos = (RAJE TBO ®pn (deg) 289.8 + 7.0 p 314.8 + 3.4 8 110.3 + 7.0 H 134.8 + 3.4 
Polarization perpendicular to the meridional plane 
u ò B (deg) 117.0 + 4.3 137.2 + 2.4 63.0 + 4.3 42.8 + 2.4 
ose = (PAJ Phn (deg) : 314.5 + 3.4 d 290.3 + 6.9 i 134.5 + 3.4 B 110.3 + 6.9 
u a B (deg) 42.8 + 2.4 63.0 + 4.3 137.3 + 2.4 117.0 + 4.3 
Coah = ERU ®pn (deg) Af 69.7 + 6.9 8 45.5 + 3.4 9 249.7 + 6.9 E 225.5 + 3.4 


measurements are not able to differentiate between inclinations torb and 180? — 25,5. The possible 
combinations result in 16 different geometrical arrangements of the black hole and orbital spins 
that satisfy the observational constraints. These 16 cases can be reduced to eight different values 
for the misalignment angle 5: four for misalignment less than 90° and four for misalignment 
between 90? and 180? for the retrograde rotation of the black hole (Table [S5). 

The probability distribution for the orbital inclination torb was assumed to be a Gaussian 
with the peak at 73° with lo error of 6° and a cutoff at 81° (72, 22). For an alternative case 
of inclination exceeding 90°, the distribution mirror reflected relative to 90° is considered. 
Other parameters are assumed to follow a Gaussian distribution with corresponding 1c errors. 
Using Monte-Carlo simulations (33), we obtain the probability distributions for 6 and ©, using 
Equations and (STT), respectively. Their mean and standard deviation are given in Table[S5] 
In Fig. ISa|we show the posterior probability distribution for 6 for the eight different cases from 
Table The probability distributions for P»n for the 16 cases from Table [S5] are shown in 


Fig. |S5| 
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Figure S4: Probability distribution function for the misalignment angle. Distributions nor- 
malized to the peak values are shown for the eight different cases presented in Table[S5| 
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Figure S5: Probability distribution function for the black hole spin azimuthal angle. Dis- 
tributions normalized to the peak values are shown for the 16 different cases (A-P) presented in 


Table 
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